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Abstract 

Recently it was shown that the Born-Infeld-type modification of the quadratic 
Yang-Mills action gives rise to classical particle-like solutions in the flat space which 
have a striking similarity with the Bartnik-McKinnon solutions known in the grav- 
ity coupled Yang-Mills theory. We show that both families are continuously related 
within the framework of the Einstein-Born-Infeld theory through interpolating se- 
quences of parameters. We also investigate an internal structure of the associated 
black holes. It is found that the Born-Infeld non-linearity leads to a drastic mod- 
ification of the black hole interior typical for the usual Yang-Mills theory. In the 
latter case a generic solution exhibits violent metric oscillations near the singular- 
ity. In the Born-Infeld case a generic interior solution is smooth, the metric has the 
standard Schwarzschild type singularity, and we did not observe internal horizons. 
Such smoothing of the 'violent' EYM singularity may be interpreted as a result of 
quantum effects. 
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Physical significance of the particle-like solutions to the Einstein- Yang-Mills (EYM) 
field equations found by Bartnik and McKinnon (BK) flj (for more details see a review 
paper ||) as well as their possible role in the string- inspired models remains rather ob- 
scure. These 'particles' where found to have a sphaleronic nature and they could be 
responsible for fermion number violating effects. However, it is not clear whether these 
purely classical solutions can survive the embedding into some quantum framework. An- 
other related puzzle is the singularity structure of the associated black holes ||, ||. It 
was shown that the metric inside the classical EYM black holes exhibits violent oscilla- 
tions near the singularity, which go far beyond the classical bounds [7|, ^|, [|. Presumably 
these oscillations must be regularized in the quantum theory, but no relevant explanation 
was suggested so far. 

To probe the effect of quantum corrections to the Yang-Mills lagrangian within the 
string/M theory framework it is tempting to utilize the Born-Infeld modification of the 



Yang-Mills action describing the low energy dynamics of D-branes JT0[. Classical solu- 
tions to both the Abelian and non-Abelian Born-Infeld theories received recently much 
attention. Although there are certain complications in the definition of the trace over 
the gauge group generators in the Born-Infeld action [[□], |12[ [13| , it is believed that the 
simplified 'square root' form of the lagrangian gives correct (at least qualitatively) pre- 
dictions concerning solitons. Magnetic monopoles were shown to persist in such a theory 
n, [T5|, IE|. A new type of a particle-like solution in the non-Abelian Born-Infeld model 



was obtained by Gal'tsov and Kerner (GK) |[7). It was shown that this theory gives rise 
to flat space classical glueballs which have a striking similarity with the BK solutions. 

More close relationship between these two particle-like configurations becomes clear 
within the Einstein-Born-Infeld (EBI) theory. As was shown recently by Wirschins, Sood 
and Kunz the GK solutions survive when gravity is added, in which case the corre- 



sponding black holes also come into play. There is a substantial difference between the 
above two families which has to be clarified, however. Both the BK and GK particles 
form the same kind of discrete sequences resulting from the 'quantization' of the param- 
eter entering the boundary conditions near the origin. But contrary to the BK case for 
which this parameter sequence is convergent to a limiting value, in the GK case the cor- 
responding sequence is divergent. So it is necessary to check more carefully whether the 
both families are continuously related indeed. 

Here we investigate the non-Abelian Einstein-Born-Infeld solitons and black holes 
in more detail. We show that within this framework one has a unique family of the 
regular particle-like solutions smoothly interpolating between the GK and BK solutions 
while the gravitational coupling constant varies from zero to a large value in units of the 
Born-Infeld 'critical field' parameter. We also investigate the interior structure of the 
non-Abelian EBIYM black holes and find that the problem of violent metric oscillations 
here is resolved indeed. 

Recall that the static spherically symmetric EYM field equations admit three smooth 
branches of local solutions near the singularity exhibiting the Schwarzschild, Reissner- 
Nordstrom and the imaginary- charge Reissner-Nordstrom type behavior. Neither of these 
three, however, has a sufficient number of free parameters to be generic. Therefore when 
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one is moving from the event horizon into the black hole interior one can meet smooth 
local solutions near the singularity at best for discrete values of the black hole mass 
0. The generic EYM black hole interior looks very differently from other explicitly 
known cases. When the singularity is approached the metric exhibits oscillations with 
an infinitely growing amplitude and an infinitely decreasing period. Clearly, the classical 
bounds are exceeded after a few oscillation cycles. In the non-Abelian EBI theory we can 
also find several local solution branches near the singularity, but the situation is drastically 
different. The local solution which has a sufficient number of parameters to be generic has 
a perfectly smooth Schwarzschild type behavior. Therefore the problem of oscillations is 
resolved in a natural way. 



As in [14, 17, 18] we assume the 'square root' form of the non-Abelian Born-Infeld 



action which, as we believe, describes well enough particle-like solutions being at the same 
time much simpler than the favored by strings symmetrized trace [1 1 action. Thus the 
EBIYM action is chosen as 

S = J (~R + 4G/3 2 (1 - K)) d 4 x, (1) 
where 



* = f + - i^i v%*r?- (2) 

Without loss of generality the dimensionless gauge coupling constant (in the units % = 
c = 1) will be set to unity, so we are left with two (dimensionfull) parameters: the BI 
'critical field' (3 of dimension L~ 2 , and the Newton constant G of dimension L 2 (Planck's 
length). From these one can form a dimensionless constant 

9 = GP, (3) 



which is the only substantial parameter of the theory ||18fl . The decoupling of gravity 
corresponds to the limit g — > 0, while the g — > oo limit (after rescaling) gives the EYM 
theory. 

We consider the case of the SU(2) gauge group assuming for the YM field a usual 
spherically symmetric static purely magnetic ansatz 

A = (1 - w (r) ) (T e sin Odtp -T^dO), (4) 

where a rotated basis for the gauge group generators is used. The spacetime metric is 
parameterized as follows: 

ds 2 = Na 2 dt 2 ~^--r 2 (d9 2 + sin 2 6d V 2 ) . (5) 

The equations of motion after a rescaling of the radial coordinate (/3) 1 ^ 2 r — ► r (making r 
dimensionless) take the form of two coupled equations for N and w: 

fNw'Y _ w(w 2 - 1) 2gw' 3 N 
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(Nr)' = 1 + 2gr 2 {\ -K), (7) 

where now 

/ m7 2 (i -w 2 ) 2 

and primes denotes the derivatives with respect to r. The third is a decoupled equation 
for a 

which can be easily solved once the YM function w is found. 

We are interested in the asymptotically flat configurations such that the local mass 
function m(r), defined through 

N=l- 2 -^, (10) 

r 

has a finite limit m — * M as r — > oo, while a — > 1. Like in the EYM case, it can be easily 
derived from the Eq. ([?]) that finiteness of M implies the following asymptotic behavior 
of the YM function: w — ±1 + (^(r" 1 ). In this limit 1Z — > 1 so the BI non-linearity is 
negligible. 

Let us first discuss the globally regular solutions which start at the origin with the 
following series expansion (its convergence in a non-zero domain may be proved by the 
standard methods |T9[|): 

6r 4 (3&(446 2 + 3) + Ag\2%b 2 + 3 - (46 2 (486 2 + 13) + 3)TZ Q 

w = l-6r 2 H * L + o(r 6 ), (11) 

30(46 2 + l) v ;> v ; 

2r 2 ^ ^ N 16^6 2 r 4 ( ^ 2 



AT = l + — g (l - Ko) - ^ - {g(n - I) 2 + SbKo) + 0(r b ), (12) 

where TZq is a limiting value of the square root at the origin 

U -> = Vl + 126 2 . (13) 

This local solution has a unique free parameter b. In the limiting case g = it coincides 
with that found in the flat space case [171], while to make contact with the corresponding 



expansion of the EYM theory one has to take the limit g — > oo with a simultaneous 
rescaling b — ► Near the origin the spacetime is flat. 

Like in the EYM case, matching of the local solution departing from the origin as 
(|IT|.|T^) and meeting the conditions of asymptotic flatness can be achieved for a discrete 
sequence of the free parameter values b = b n (g) labeled by the number n of zeroes of w(r). 



The proof of existence may be given along the lines of JT9], [L7|] , here we do not enter into 
mathematical details concentrating rather on the qualitative physical picture. 

We have investigated numerical solutions in a large range of g. Typical behavior for 
small g is shown of Figs. (|1],0). The YM curves for g = .001 are practically undistin- 
guishable from those found previously by Kerner and one of the authors [0 in the flat 



space BIYM model. The region of oscillations corresponds to an unscreened Coulomb 
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charge inside the particles. One can see that this region expands both in the direction 
of small and large r with growing n. The metric deviate from the flat one rather weakly 
for small g (weak gravity). One observes a stabilization of the metric in the region of the 
w-oscillations. 

The amplitude of oscillations of w in the intermediate region decreases, so in this region 
the solution can be regarded as approaching the embedded Abelian solution w = 0. As 
was observed by Wirschins, Sood and Kunz [n|, the metric for small g also approaches 
the corresponding Abelian Bion metric |2(J (with zero 'seed' mass in the singularity). To 
avoid confusion it is worth noting that for odd n the Yang-Mills topology of the EBIYM 
regular solutions (kink) is essentially different form that of the embedded Abelian one 
w = (trivial). Thus it would be misleading to say that the sequence of non- Abelian 
EBIYM solutions converges to an Abelian one in the global sense. 

For large g the regular EBIYM solutions after the coordinate rescaling r — > ^fgr tend 
to the BK solutions of the EYM system, see Fig. [3] for (7 = 10. Comparing this with the 
weak gravity behavior (Fig. H) we observe that gravity reduces the region of the unscreened 
charge. The metric deviates substantially form the flat one. The function N(r) for large 
n exhibits a deep well (an 'almost' horizon), but remains always strictly positive. 

The sequences b n of discrete values of the parameter in the expansion (|TTJ) behaves 
very differently for small and large g. As was found in JL7], b n in the flat space are quite 



big with respect to the corresponding BK values, and the sequence does not converge 
with growing n. Contrary to this, the BK sequence b n rapidly converges to a limiting 
value &oo. In this case there is an additional limiting solution with different space-time 
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We have obtained numerically the interpolating functions b n (g) for several n clearly 
demonstrating that these two extrema are continuously related indeed (Fig. |5|). For small 
values of g these functions tend to the GK values |L7| 



MO) = b™, (14) 

while in the opposite limit g — > 00 one recovers the convergent sequence of rescaled BK 
values H, 1 

b* K = ]imb n (g)/g. (15) 

y * CXJ 

The corresponding solutions of the EYM system with the standard YM lagrangian are 
recovered in terms of the rescaled coordinate r = r / yj~g. In the intermediate region all 
parameter functions b n (g) were found to be monotonously varying between the above 
extrema. 

Masses of the regular solutions as functions of the effective gravitational coupling are 
shown on Fig. || For vanishing g one recovers the GK masses after a rescaling 

\imM n {g)^M^ K g. (16) 
We recall [[17]] that the sequence M^ K converges to the mass of an embedded Abelian 



solution. The reason is that the main contribution to the mass comes from the region of 
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oscillations where w approaches an Abelian value w = with growing n. But for any n 
the limiting values w(0),w(oo) are equal to ±1, so, as we have already noted, the global 
topology of no n- Abelian solutions is entirely different. In the limit of strong gravity one 
recovers the BK masses after a rescaling 



Jim M n {g) - M™y/g. 



(17) 



Now discuss the black holes. These are parameterized by the horizon radius and 
the value of the YM function Wh = w(rh) at the horizon. The series expansions near the 
regular horizon reads: 



w 
N 



K 



^^(r-r h ) + 0((r-r h r) 



Wk+ r h N h 
N' h (r-r h ) + 0((r-r h ) 2 
1 



18) 



1 + 2g\ 2 h 1 + Jl + (wl - iy/ri 



Asymptotically flat solutions with such boundary condition are likely to exist for any hori- 
zon radii r^. Exterior black hole solutions are very similar to regular solutions, especially 
for small r^. So our main interest is in the interior solutions. We start by listing various 
series solutions that can be obtained near the singularity. 

For 



Let us first explore the series expansion for an embedded Abelian solution [20, 21 



the unit magnetic charge (what corresponds to w = 0) the local mass function satisfies 
the equation 

g [Vr 



m 



V 4 + 1 



(19) 



Expanding the square root at small r and integrating one obtains: 



2m n 2 

N = ° + 1 - 2^7 + -gr< 

r 6 



1 4 

-or 
5 y 



Ofr 5 



(20) 



The 'seed' mass parameter mo may be positive, negative or zero. Positive m corresponds 
to a timelike singularity of the Schwarzschild type. Negative m corresponds to a space- 
like singularity, in which case an internal horizon also exist (though contrary to a more 
common example of the Reissner-Nordstrom metric, the 'local charge' term Q 2 /r 2 now 
is absent). Vanishing mo is particularly interesting. The metric near the origin is then 
locally flat unless g = 1/2 in which case the limiting value N(0) shrinks to zero. For other 
values of g there is a conical singularity PD| , the critical value g = 1/2 corresponding 
to an extremal deficit angle. Within the present framework the embedded Abelian black 
holes correspond to an identically vanishing function w. It can be shown that the local 
series solution starting at r = with zero initial value Wq = and arbitrary m generates 
this global embedded Abelian solution. 

For non-Abelian solutions the value of the YM function at the singularity w(0) = Wq 
should therefore be non-zero. We have found the generalization of the series expansion 
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(pO|) with wq ^ 0. It is valid for the non-zero seed mass parameter itiq: 



W = W " 1+ 2^ + l«l 3(2e - 1) - 49c| +0(r) ' (21) 



2m 3qew?,r 9 / 2 eiOn , . i 

eJV = + 1 - 2 S e + -|-JL 4- 9 r* - + ^ [3(5e - 2) - S S e] ) + O(r'), 



where e = 1 — iOg. This local solution has two free parameters w ,m . For u> = it 
coincides with (p0|). If u> 7^ becomes singular in the limit m = 0. The search for local 
solutions with tuq = shows that in this case either wq = in which case we come back to 
the Abelian embedded solution w = 0, or Wq — ±1, then we recover the regular solution 
(0.0). It is also worth noting that, contrary to the EYM case, in the EBIYM theory 
there are no local solutions with iV ~ r~ 2 behavior at the singularity (Reissner-Nordstrom 
type). 

Since the asymptotic solution fl2"T|) fails to contain a sufficient number of free param- 
eters to be a generic solution, (this number is equal to three for the system of equations 
(0,[7|)), the question is how a generic solution looks like near the singularity, in particular, 
whether it admits any series expansion. In the case of the ordinary Yang-Mills lagrangian 
such an expandable solution does not exist at all, one finds that a generic solution has a 
non-analytic oscillating behavior J?). Here the situation is different, although the generic 
local solution around the singularity still exhibits non-analyticity in terms of the variable 
r. It turns out to be series expandable but in terms of the y/r: 

r- a 2 w r a (24g 2 a 2 - 32g 2 a 2 e - 3c 2 ) r 3 / 2 
w = w + aV-r- — + 

a 2 ((16o 2 a 2 — 32o 2 a 2 e — 15c 2 ) w — 8q 2 cae) r 2 -/,. 

— 32?^ — + 0{r 

N = 1 - 2^- + cV~r - " (3l "° C + " g2e) r + Q(r 3 ' 2 )- (22) 
err e 

This local solution contains three free parameters wq, a, c. The singularity is of the 
Schwarzschild type, the seed mass mo is strictly positive. Two leading terms in expansion 
of 'kinetic' (negative for N < 0) and 'potential' (positive) terms in 7Z cancel, so the leading 
behavior in singularity is 

4^3/2 

We have tested for various rh,g,n that a continuation of the exterior black hole solu- 
tions under the horizon meets this generic asymptotic solution indeed. In all numerical 
experiments the function TV remained negative under the horizon and no internal Cauchy 
horizons were met. Typical global black hole solutions are shown on the Figs. |7||||9| for 
Th = 1, g = 1, n = 1, 2. The YM function w outside the horizon has qualitatively the 
same behavior as in the regular case. Inside the horizon it remains perfectly smooth and 
tends to a finite limit w ^ at the singularity. Recall that for the ordinary quadratic 
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YM lagrangian the function w inside the horizon of the EYM black holes has a rather 
sophysticated behavior: while w itself tends to a finite limit w as well, its derivative 
exhibits a sequence of infinitely increasing absolute values at tiny intervals, whose length 
tends to zero |7], Such a behavior causes oscillations of the local mass m(r) with an 
infinitely increasing amplitude. At the beginning of each oscillation cycle the metric func- 
tion N(r) takes values very close to zero ('almost' Cauchy horizons), then an exponential 
growth of m(r) starts. After a few oscillation cycles the maximal values of m attained in 
subsequent cycles become of the googolplexus order, obviously lying beyond any classical 
bounds. Contrary to this, in the EBIYM case we observe a smooth m(r) inside the hori- 
zon up to the singularity where m(r) has a finite positive value (Fig. ||). Similarly, the 
second metric function a tends smoothly to a finite value at the singularity (Fig. |9|) 

We conclude with the following remarks. The BK solutions were found for the gravity 
coupled Yang-Mills theory with the usual quadratic lagrangian. Now they were shown to 
be a strong gravity limit of the gravity coupled BIYM theory which can be interpreted 
as an effective YM theory including string quantum corrections. Within this theory there 
exists a limit of decoupled gravity, in which qualitatively similar solutions continue to 
exist. This shows the way how the BK solutions can be incorporated into the (quantum) 
string theory. 

A particularly interesting implication of this reasoning is the resolution of the prob- 
lem of 'violent' oscillating singularities typical for hairy EYM black holes. Born-Infeld 
corrections perfectly regularize the behavior of the metric near the singularity. It is worth 
noting that the generic singularity is timelike, in conformity with the strong cosmic cen- 
sorship hypothesis. In numerical experiments we did not observe internal horizons. In 
principle, such horizons could emerge in pairs, this question is worth to be investigated 
in more detail. 

This work was supported in part by the RFBR grant 00-02-16306. 
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Figure 1: First five solutions w n (r) for weak gravity (g = .001). The oscillation region, 
which corresponds to the localization of an unscreened magnetic charge inside the BIYM 
particle, expands with growing n. The amplitude of the first and the last oscillation cycles 
remains roughly the same, while the amplitude of the intermediate cycles tends to zero 
with increasing n. All these curves practically coincide with those found in the flat space 
BIYM theory by Gal'tsov and Kerner. 
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Figure 2: Metric function N(r) for regular solutions with g = .001. For all solutions the 
deviation of the metric from a flat one is within 2%. For higher node numbers n one 
observes a stabilization of iV near the minimal value, this corresponds to the oscillation 
region of w(r). 
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Figure 5: Dependence of the parameters b n for regular solutions on the effective gravity 
coupling constant g (in logarithmic coordinates). The left side of the picture correspond 
to the GK limit. With growing g higher-n curves merge, this is reminiscent of the rapid 
convergence of the sequence b„ K = lim g ^ 00 b n (g)/g. 
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Figure 6: Dependence of masses of regular solutions on the gravity coupling g 
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Figure 7: YM function w(r) for n = 1, 2 EBIYM black holes with r h = 1.0, # = 1.0,. At 
the singularity these curves tend to constant values non-equal to unity. 
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Figure 8: The mass function m(r) function for the EBIYM black holes with — 1.0, (7 = 
1.0, n = 1,2. Inside the horizon m(r) remains a smooth function tending to a constant 
limit at the singularity. There are no internal horizons, the singularity is timelike. 
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Figure 9: The metric function a{r) for the EBIYM black holes with r h = 1.0, g = 1.0, n = 
1,2. 
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